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Abstract 

We have searched for the C-violating decay n — > 777 in a sample of ~18 million 
77 mesons produced in <f> — > 777 decays, collected with the KLOE detector at the 
Frascati ^-factory DA$NE. No signal is observed and we obtain the upper limit 
BR(r/ -> 777) < 1.6 x 10~ 5 at 90% CL. 



The decay 77—^777 is forbidden by charge-conjugation invariance, if the weak 
interaction is ignored. The present limit for the 77 — > 37 branching ratio, 
BR(rj -> 37) < 5 x 1(T 4 at 95% CL, is based on the result of the GAMS2000 
experiment at Serpukhov [1], which studied neutral decays of 77 mesons from 
the reaction ir~p — > rjn at a beam momentum of 30 GeV/c. 

We have searched with KLOE for the decay 77—^777 among four-photon events, 
corresponding to the two step process <fi— ^777, 77^777. The KLOE detector 
[2-5], operates at the Frascati e + e~ collider DA$NE [6], which runs at a 
CM energy W equal to the 0-meson mass, VF~1019.5 MeV. Copious 77-meson 
production is available from the decay — > 777, with a branching ratio of 1.3%. 
The highest ^-production rate that has been obtained to date was ~240 0/s, 
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corresponding to ~3.1 77/s, in October 2002. At DA$NE, because of the beam- 
crossing angle, mesons are produced with a small transverse momentum, 
12.5 MeV/c, in the horizontal plane. The present analysis is based on data 
collected in the years 2001 and 2002 for an integrated luminosity of 410 pb" 1 , 
corresponding to 1.8 x 10 r\ mesons produced. 

The KLOE detector consists of a large cylindrical drift chamber [2], DC, sur- 
rounded by a lead/scintillating- fiber sampling calorimeter [3], EMC, both im- 
mersed in a solenoidal magnetic field of 0.52 T with the axis parallel to the 
beams. Two small calorimeters [4] are wrapped around the quadrupoles of 
the low-/3 insertion to complete the detector hermeticity. The DC tracking 
volume extends from 28.5 to 190.5 cm in radius and is 330 cm long, centered 
around the interaction point. The DC momentum resolution for charged par- 
ticles is 5p±/p±=0A%. Vertices are reconstructed with an accuracy of 3 mm. 
The calorimeter is divided into a barrel and two endcaps, and covers 98% of 
the total solid angle. Photon energies and arrival times are measured with 
resolutions a E /E = 0.057/y^ (GeV) and a t = 54 ps/^/E (GeV) © 50 ps, 
respectively. Photon-shower centroid positions are measured with an accuracy 
of a — 1 cm/ J E (GeV) along the fibers, and 1 cm in the transverse direc- 
tion. A photon is defined as a cluster of energy deposits in the calorimeter 
elements that is not associated to a charged particle. We require the distance 
between the cluster centroid and the nearest entry point of extrapolated tracks 
be greater than 3xa(z, <fi). 

The trigger [5] uses information from both the calorimeter and the drift cham- 
ber. The EMC trigger requires two local energy deposits above threshold 
(E > 50 MeV in the barrel, E > 150 MeV in the endcaps). Recognition and 
rejection of cosmic-ray events is also performed at the trigger level by check- 
ing for the presence of two energy deposits above 30 MeV in the outermost 
calorimeter planes. The DC trigger is based on the multiplicity and topology 
of the hits in the drift cells. The trigger has a large time spread with respect 
to the time distance between consecutive beam crossings. It is however syn- 
chronized with the machine radio frequency divided by four, T sync =10.85 ns, 
with an accuracy of 50 ps. For the 2001-2002 data taking, the bunch crossing 
period was T=5.43 ns. The time (To) of the bunch crossing producing an event 
is determined offline during event reconstruction. 

The sensitivity of the search for r/^777 in KLOE is largely dominated by 
the ability to reject background. The dominant process producing four pho- 
tons is e + e~ — > w-f, due to initial-state radiation of a hard photon, followed 
by uj — > 7r°(^ 27)7. Other processes with neutral secondaries only are also 
relevant. They can mimic four-photon events because of the loss of photons, 
addition of photons from machine background, or photon shower splitting. All 
the above effects are very difficult to reproduce accurately with Monte Carlo 
(MC) simulation. We therefore base our background estimates on data, and 



3 



use the MC only to evaluate the efficiency. An 77— >37 generator using phase 
space for the internal variable distribution in the three-body decay has been 
used to produce 120,000 — > 777, 77 — > 37 events. 

For the analysis, only events without charged particle tracks are considered. 
The central value of the position of the beam-interaction point (IP), the CM 
energy, and the transverse momentum of the are obtained run by run from 
large samples of Bhabha scattering events. The following requirements have 
been used to isolate 0-^47 candidates: 

(1) The four photons must have 

• reconstructed velocity consistent with the speed of light, \t — r/c\ < 5cr t , 
where r is the distance traveled, t is the time of flight and a t is the time 
resolution; 

• photon energy E 1 >50 MeV; 

• photon polar angle 9 >24.5°. 

(2) The total energy and momentum of the four prompt photons must satisfy 
E i £7 i >800 MeV and |E^| <200 MeV/c; 

(3) The opening angle between any photon pair must satisfy 77 > 15°. 

83,906 events pass the cuts above. A kinematic fit is used to improve the 
energy-momentum resolution. The input variables Xi of the fit are 

- the coordinates of the photon clusters in the calorimeter; 

- the energies of the clusters; 

- the times of flight of the photons; 

- the coordinates of the e + e~ interaction point; 

- the energy and momentum of the meson. 

We minimize the y 2 function 

i a i j 

where Fj(fik) are the energy, momentum, and time constraints and A,; are 
Lagrangian multipliers. The x 2 value of the fit is used to reject background. 
Events with \ 2 < 25 are retained, the number of degrees of freedom being 8. 

After this cut we are left with 52577 events. The residual background after 
the cut is due to events with neutral pions (fig. 1, left), coming mainly from 
e + e~ — > with uo — > 7r°7. This can be seen in fig. 1, right, where the invariant 
mass of the tt° and the highest energy photon, in the 7r°77 hypothesis, 7hi, 
shows a clear peak at the u> mass. Other background sources with a 7r° in the 
final state are the decays 0— >7r°7, — > f j — > 7r°7r°7, and — > a 7 — > r]7r 0r y. 
We reject the main part of these events by a cut on the invariant mass of any 
photon pairs: 90 < 774(77) < 180 MeV. 
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Fig. 1. Distribution of the invariant masses 771(77) computed for all photon pairs, 
left and m(7r°7hi), right. The shaded interval is removed before further analysis. 

8,268 events survive the cuts. In the decay 0^r/7, the energy of the recoil 
photon in the CM of the is 363 MeV. In the complete chain 0^7/7, 7/^37, 
363 MeV is also the most probable energy of the most energetic photon, 7^. 
Fig. 2, left, shows an MC simulation of the 7hi energy spectrum for the signal. 
Fig. 2, right, shows the £"(7^) distribution for the data sample. No peak 
is observed around 363 MeV. To evaluate an upper limit on the number of 




Fig. 2. Distribution of the energy E^^i), in the 4> CM, for the MC simulated signal 
(left) and for the data (right). See text for discussion of the background fit. 

77 — > 37 events, we choose as the signal region the interval 350< -E^m) <379.75 
MeV (17 bins, 1.75-MeV wide). We estimate the background by fitting the 
£"(7hi) distribution on both sides of the expected signal region, in the intervals 
280< E <350 and 379. 75< E <481.25 MeV. We fit the background using 3 rd 
to 6 th order polynomials. The 5 th order polynomial shown in fig. 2 gives the 
best fit, with x 2 /dof = 78/92 = 0.85. We use the result to obtain the expected 
number of background events in each bin, Nf. The total number of observed 
events in the signal window is 1513 while from integration of the polynomial 
we obtain 1518 events in the same region. 

The upper limits have been evaluated using Neyman's construction procedure 
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[7]. To evaluate the agreement with the background distribution in the signal 
region, we use 

F _ y (^-iVf) 2 
r 2^ N b 

i i 

where AT* is the number of observed counts in the i th bin, and the sum is over 
bins in the signal region. We obtain the distribution function for F for various 
values of the number of signal counts s as follows. First, we construct the values 
Ni by sampling a Poisson distribution with mean (iVj(s)) = N^ + s x fa, where 
fi is the fraction of signal events (fig. 2, left) in the i th bin, and evaluate F. 
Repeating this procedure 10 6 times for each value of s then gives the complete 
p.d.f., which is numerically integrated to obtain the 90% and 95% contours in 
Neyman's construction. We then evaluate F using the observed N{. We find 
F = 13.45, from which we obtain 

< 63- 1 at 90% CL ; < 80 - 8 at 95% CL - 

To convert this result into an upper limit for the branching ratio, we nor- 
malize to the number of rj — > 3ir° events [8] found in the same data sample, 
N(r] -> 3tt°)=2,431,917. The efficiencies are e(r) -> 3vr )=0.378±0.008(s?/st.)± 
0.001(stat.) and e{r) -> 37) = 0.200 ± 0.001(stat.) ± 0.002(syst.) ± 0.006(xL)- 
The systematic error includes residual uncertainities on the photon detection 
efficiency [9]. For the ratio of the two branching ratios we obtain 

*4| = N <r* f "- 3 "' < 4.9 x 10- 5 90% CL, 
- 3-°) e „_ 37 " 6 3 x iQ _ 5 m% cl 

Using the value BR(r/ -> 3tt°) = (32.51 ± 0.29)% [10], we derive the upper 
limit 

BR(r] -> 37) < 1.6 x 10~ 5 at 90% CL and < 2.0 x 10~ 5 at 95% CL. 



The efficiency quoted above for 77— >37, which depends on the cut x 2 < 25 
applied after the kinematic fit of all four-photon events, is evaluated by MC 
simulation. We check the validity of the MC result by comparison with the 
X 2 distribution for radiative events e + e~— ^uo— >777r°— ^47. A sample of these 
events is selected from among all four-photon candidates by requiring 128 < 
771(77) < 145 MeV for the neutral pion and 760 < m(7r°7) < 815 MeV for 
the u. The fraction of these events with x 2 < 25 after the kinematic fit differs 
from the MC estimate by ~3%. This value is included in the quoted error for 
e(rj -> 37). 

To check whether the kinematic fit introduces a bias in the energy distribu- 
tion of the signal photons, we have analyzed a sample of (ft— >t]'~f TCC ^777 rec 
events, in which the energy of the recoil photon is the same as in the case of 
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interest. Fig. 3 shows the energy distribution of the photons as obtained after 
the kinematic fit for data and MC events. The two distributions are in good 
agreement within errors. 



4000 
3000 
2000 
1000 



" 350 355 360 365 370 375 
E(Y) MeV 

Fig. 3. Distribution of E{^) for the (ft — > 777 — > (77)7 event sample for the data 
(points) and MC simulated events (continuous histogram). 

The stability of the upper limit versus the background estimate has been 
checked by comparing the results of polynomials of different degree for fit- 
ting the E("f\a) distribution outside the signal region. A 3 rd order polynomial 
doesn't describe the background shape well. A 4 th order polynomial gives a 
lower value for the signal yield, while a 6 th order polynomial gives the same 
result. We have also checked the stability of the result by changing the win- 
dow chosen for evaluation of the upper limit obtaining a maximum variation 
of 11%. We have also evaluated the 77 — > 37 acceptance using the matrix 
element of ref.[12] and we find a value 5% lower. Therefore systematic ef- 
fects can be summarized: background estimation and window variation 11%, 
e(r) — > 37r°)/e(?7 — > 37) 1%, \ 2 cu f 3%, decay model 5%. We thus feel confident 
about the procedure adopted. Our limit 

BR(?7 -> 777) < 1.6 x 1CT 5 at 90% CL or < 2.0 x 10~ 5 at 95% CL 

is the strongest limit at present against possible violation of charge-conjugation 
invariance in the decay 77^37. 2 An estimate for F(r) — > 37), including con- 
tributions from weak interactions, is given in Ref. 11. Using the estimate for 
7T° — > 37 [12], one finds BR(?7 — > 37) < 10~ 12 , which is quite a long way 
from the result above. The absence of the decay 77— >37 therefore confirms 
the validity of charge-conjugation invariance in strong and electromagnetic 
interactions. 

We thank the DA$NE team for their efforts in maintaining low background 

2 A preliminary unpublished 2002 result by the Crystal Ball, BR(r] —> 37) < 1.8 x 
10" 5 at 90% CL, is mentioned in [13]. 
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